Quantum Antiferromagnetism of Fermions in Optical Lattices with Half-filled p-band 
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We study Fermi gases in a three-dimensional cubic optical lattice with five fermions per site, i.e. 
the s-band is completely filled and the p-band with three-fold degeneracy is half filled. We show that, 
for repulsive interaction between fermions, the system will exhibit spin-3/2 antiferromagnetic order 
at low temperature. This conclusion is obtained both in strong interaction regime by strong coupling 
expansion, and in weak interaction regime by Hatree-Fock mean-field theory with analysis of the 
Fermi surface nesting. We also show that, in the strongly correlated regime, the Neel temperature 
for p-band antiferromagnetism is two to three orders of magnitudes higher than that of s-band, 
which is much more promising to be attained in cold atom experiments. 



I 



I 

o 



> 

oo 
m 



o 



I. INTRODUCTION 

Studying strong correlations in lattice Fermi gases has 
now become an emerging forefront of cold atom physics, 
mainly because of the flexibility of varying interactions 
and controlling filling number to access various quantum 
phases. One major effort pursued in many labs now is 
to use ultracold Fermi gases to simulate Hubbard model, 
which exhibits various metallic, insulating and supercon- 
ducting phases, with different magnetic behaviors^. De- 
spite of huge efforts extending over several decades, there 
are still lots of controversial issues on this model, which 
hopefully cold atom experiments can shed light on. 

The most solid conclusion of Hubbard model is drawn 
for a half- filled cubic lattice which says a spin- 1/2 anti- 
ferromagnetic (AF) order exists for all range of repulsive 
interaction. For weak interaction, the AF order arises 
from the nesting geometry of Fermi surface, while for 
strong interaction, it is caused by the super-exchange 
interaction^ii^i^i^. Observing this AF order in cold atom 
experiments, especially in strong interaction regime, is 
considered as a hallmark of seeing strong correlations in 
lattice Fermi gases and the starting point for more am- 
bitious goals along this direction. 

On the other hand, in optical lattices physics can be 
much more richer than the single-band Hubbard model. 
When there are more than two fermions at one lattice 
site, fermions will start to occupy excited bands like the 
p-band. Unlike in condensed matter systems, there is 
no Jahn- Teller effect in optical lattices, so the three-fold 
orbital degeneracy of the p-band is well maintained. 

Recently there are increasing theoretical and experi- 
mental interests in studying bosons in the p-band, where 
a host of intriguing quantum phenomena due to orbital 
degeneracy has been pointed oulj^ii^iiSiiiS. We expect the 
physics of p-band fermions will be more easier for ex- 
perimentally study, because the Pauli exclusion princi- 
ple prevents the decay into the lower band. Thanks to 
the orbital degeneracy, p-band Fermi gases are expected 
to exhibit more diverse phases comparing to single-band 
Hubbard model. Nevertheless, the first step toward re- 
vealing these exciting phases is to understand the un- 



frustrated spin order at half-filling. Here we consider a 
three-dimensional cubic lattice whose three degenerate 
p-bands are half filled by three fermions per site. (The s- 
band is already completely filled by another two fermions 
and do not need to be worried about.) 

With the complexity in p-band model, even for half- 
filling there are issues not clear, for instance (i) whether 
the p-band fermi-surface still has some nesting properties 
which give rise to a magnetic instability in weak interac- 
tion regime; (ii) with multiple hopping channels in differ- 
ent orbitals and multiple choices of intermediate states, 
whether super-exchange processes can still yield a simple 
Heisenberg-type model as in the s-band situation; and 
(iii) as there are various types of orders due to orbital 
degrees of freedom, whether the strong and weak inter- 
action regime share the same order or there is quantum 
phase transition in between. These questions will be an- 
swered in this work. The main results of this paper can 
be summarized as: 

(I) In strong interaction regime, at each isolated site 
the ground state of three fermions in p-orbitals are spin- 
3/2 states with four- fold degeneracy due to the Hund's 
rule. The coupling between two neighboring spins is 
caused by virtual hopping of fermions, which gives rise 
to an isotropic spin-3/2 Heisenberg model: 



H,7 — ^cx^ 



SiSj —Jiinj 



(1) 



where Si is an on-site spin-3/2 operator. Jox > which 
means the ground state of this system has spin-3/2 anti- 
ferromagnetic order. 

(II) In weak interaction regime, the Fermi surface of 
a half-filled p-band exhibits a perfect nesting symmetry, 
with the nesting momentum tt — (tt, tt, tt). Mean- field 
analysis shows this nesting will induce antiferromagnetic 
order at low-temperature, and the order parameter coin- 
cides with the 5 = 3/2 spin order in strong interaction 
regime. We also show that this order should be the only 
order parameter. 

Based on (I) and (II), we anticipate that, despite the 
complexity with p-band orbital degeneracy, the half-filled 
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p-band Hubbard model exhibits only the spin-3/2 antifer- 
romagnetic order for all range of interaction strength, and 
the system undergoes a crossover from weak to strong in- 
teraction, a situation similar to the half-filled single-band 
Hubbard model. 

In the experiments of cold atoms in optical lattices, the 
major difficulty in achieving AF order in strongly corre- 
lated regime comes from the fact that the Neel temper- 
ature is very low because of the ultra-small energy scale 
of super-exchange interaction. In this work we show that 
in strong interaction regime the Neel temperature for p- 
band AF can reach a few nano-Kelvin, which is two to 
three orders of magnitude higher than that of s-band. 
Because of this we suggest to study p-band Fermi gas as 
a more realistic and promising route to reveal strongly 
correlated physics in optical lattices. 

II. MODEL 

The model for Fermi gases in optical lattices contains 
a nearest neighbor hopping Ht and an on-site interaction 
H^/. The hopping term Ht can be written as: 

a,/3— a;,y ,2: 

Ht= E ^r4.,.^.+/5,p„,. + h-c., (2) 

Since optical lattice is a separable potential in three di- 
mensions (Vo(sin^ Kx + siv? Ky + s\v? Kz j), the hopping 
term conserves orbital index due to the orthogonality of 
wannier wave functions, and the hopping amplitude is 
anisotropic Tg° = tdap — ^'(1 — Sap) with t,t' > 0. And 
there is also no next nestest neighbor hopping*. Fig. 
[T]JA) is a schematic of hoppings in the p-band. t is much 
larger than t' , for instance, for an optical lattice with the 
depth Vo = 11£^R (£^r = h?K'^/{2m) is the photon re- 
coil energy), t « 20i'. For ultra-cold neutral gases, the 




(A) (i) K 



interaction is dominated by short-range s-wave repulsion 
^(ri2) — {^■nh^as/m)5{rx —^2), which takes place be- 
tween unlike spins only^. The on-site interaction Hj/ 
then has the form 

U = ^ / d^r^ r 4 4 

m J 

W^^^ fdh\w,Ar)\'\w,Ar)\' (5) 

Til J 

where Wp^ is the pa-orbital wannier wave function. There 
are totally four terms due to s-wave interaction. The first 
term is the interaction within the same orbital, the sec- 
ond term is the interaction between unlike orbitals, the 
third term represents the spin exchange between unlike 
orbitals and the fourth term represents two fermions hop- 
ping processes from one orbital to another. Usually U is 
three times larger than W. 

III. STRONG INTERACTION REGIME 

When the on-site interaction is dominative, we shall 
first ignore the hopping term and calculate all the eigen- 
states of an isolated site. The results for three, four and 
two fermions per site are summarized in following Table 
HI [Til and mil As one can see from Table HI the low-energy 
Hilbert space is constituted by four spin-3/2 states. Be- 
cause these spin states are symmetric superpositions of 
three atoms, their spatial wave functions have to be anti- 
symmetric and the interactions are completely quenched. 
In this limit, the system is a spin-3/2 Mott insulator— . 



E^O{S = 3/2) 


E = 3W{S^ 1/2) 


E>U 


1 TTT) 
1 iii) 

(1 iTT) + 1 Tit) + 1 TTD/v^ 
(UiT) + ini) + UTi)/v^ 


(1 TTi)-l TiT))/v^ 

(1 m-\tim/V2 

(1 iiT)-UTi))/^ 

(1 ni}-\iU})/V2 


double 
occupied 
states 



TABLE L Eigenstates and eigenenergies of three fermions at 
an isolated site. In these three tables, we define that, for 
example, |d ^ 0) means is double-occupied, py is occupied 
by a spin-t fermion, and pz is empty. 



FIG. 1: (color on-line) (A) Schematic of fermion hopping in 
p-band. For a fermion in pa, orbital in i-site, its hopping 
amplitude to i+x site is t, to i+y-site and i+z site is t'{t' <^ t), 
and it can not hop to Py (or p^) orbital of its neighboring site. 
(B) Illustration of the p-band Fermi surface nesting at half- 
filling. Red line and blue line are Fermi surfaces in — ky 
plane for kz = n/3 and kz — ~2tt/3, respectively. These 
two Fermi surfaces parallel with the nesting momentum 7? — 
(tt, 7r, 7r)(the solid black arrow). 



The coupling of spins between neighboring sites can be 
obtained by taking the hopping term as the second-order 
perturbation with the Brillioun-Wigner approximatiort^: 

Hj = - V PalltPe^VentPG, (6) 

where means the projection into excited Hilbert space, 
and Vg means the projection into S = 3/2 Hilbert space. 
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E = U + 2WiS =1) 


E = U + AW{S = 0) 


MTT), 1 Tdt), 1 TTd), |dii), 

\idl),\ lid), (|dn) + |diT))/V2 

(1 Trfi) + UdT))/^ 

(1 Ti d) + 1 it d))/V2 


(|dn)-MiT))/%/2 
(1 Tdi>-UdT))/V2 
(1 n d) - 1 it d))/V2 



TABLE II: Eigenstates and eigenenergies of four fermions at 
an isolated site. 



£ = 0(5=1) 


E = 2W{S = Q) 


E>U -W 


|OTT), ITOT), ITT o), 
10 ii), lioi), Uio) 
(|0Ti) + |0iT))/V2 
(1 T i) + Ii T))/v^ 
(ITT 0) + |TT0))/V2 


(|0TT)-|0TT))/V2 
(1 T T) - 1 T T))/V2 
(1 TT 0) - 1 TT 0))/^/2 


double 
occupied 
states 



TABLE III: Eigenstates and eigenenergies of two fermions at 
an isolated site. 



To proceed we make following observations: (i) The vir- 
tual hopping processes only take place within the same 
orbital, otherwise the final state will have double occu- 
pancy, which should be projected out. For instance, if a 
Px fermion hops from i-site to j-site, it must be followed 
by a Px fermion hops from j-site to i-site. (ii) By using 
Table HI [Til and Hill one can show that only one type of ex- 
cited states can be connected to the unperturbed ground 
states via hopping, which is the lowest energy states of 
four fermions (S* = 1 and E = U+2W) in one site and the 
lowest energy states of two fermions {S = 1 and E — 0) 
in its neighborhood^^. So the intermediate states always 
has energy U + 2W. 

Because of (i-ii), Hj on the link between i and i + x 
site reads 



E 



U + 2W 



E4p, =C,,„„.vC, 



Note that 



3 I E^P^.s*^* 



,Pl3 



Pf3 



where 



and with the help of Table HI one can show that 



'Pg^^Lp^'Pg = 



'i,Pl3,s' ■ 



(7) 



(8) 
(9) 



Then Hj on this link can be written as 

2|r|'°|2 

9(U + 2W) 2^ Z^^^'Pfi'^^^^'Pi^'''' 2- 

Pa SS' Pf, pjj 



)q ss' Pfj 



Pa ,P/3 S,S' 
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where 

Jox = 4(i2 + 2i'2)/(9f/ + 18W) 
We now introduce Si as 

^ 9 E E 4p, 



^(7 ss' ^i,Pa ,s' 



(10) 



(11) 



Pa ss' 



and a are the Pauli matrices. Restricted in the low- 
energy Hilbert space. Si acts as S = 3/2 spin operator. 
Hence, we have shown the first main result that 

In the strong interaction regime, the super-exchange pro- 
cesses induce an isotropic spin-3/2 Heisenberg model 



Hj 



SiSj 



1 



-niUi 



(12) 



which describes the low-energy spin dynamics of p-band 
Mott insulator. Jex > 0. Below Neel temperature, it 
will give rise to an antiferromangetic long range order of 

(5+) = (-l)i(5+)|^0 



IV. WEAK INTERACTION REGIME 

In this regime we use usual Hartree-Fock mean-field 
theory, for which the Hamiltonian is rewritten into mo- 
mentum space. The hopping term can be written as: 



(13) 



where 



^k" ~ 2tcoska — 2t' cos/c/j — /i. (14) 

For half-filling, /i = 0, and the Fermi surface for one of 
the p-bands is plotted in Fig. [HB). 

For performing mean-field decomposition, we recon- 
struct the interaction term H[/ as -I- H'^ , where 



k,k',q 



(15) 



ii'u - E (-y^^i^^"^'^ 2w±'^^±'^Q (16) 

k,k',g 

with Vl = 3(C/ + 2W) and V2 = 2{U - W)/3. Here we 
introduce 



^k,q — E ^Pa,Pa,k/'^' 



(17) 



P/3 



4 



and 



pa/3t 



,k '5'p^tp^,k)/2, 



(18) 
(19) 



4,p„,T'^k+q,p3,i is the spin lifting op- 



where S^^ ^^ ^^ 
erator. As illustrated in Fig. IHB), each of three 
p-band Fermi-surfaces exhibits nesting symmetry with 
the nesting momentum tt = (7r,7r,7r) at half-filling. It 
means that for each k at Fermi surface, k -I- tt also lo- 
cates at Fermi surface, therefore, the maximum contri- 
bution to Hf/ comes from q = tt. For this reason we 
only consider q = 7? component in the mean-field ap- 
proximation, and we introduce three order parameters 
which are A = l^i Ek(Ak..), F = 1^2 Ek(rk!l>/3 and 

E = 2iy Ek(^k'7r)' where F and S are independent of 
index a, (3 because of the permutation symmetry between 
orbitals. Note that A is the same spin order as in strong 
interaction regime, we separate the mean-field Hamilto- 
nian into two parts as Hmf = Hq -f H' where 



Ho = E 



E^k°4,p„,.Ck,p„,. + (AAj^-Hh.c.) 



A2 



and 



H' = E [s^f ^ + h.l 



V2 ^ 2W 



Note that H' does not contain f-terms because they all 
cancel out due to the permutation symmetry of three 
orbitals, therefore F = 0. 

Because the three orbitals are now separable in Hq, it 
can be easily diagnolized by Bogoliubov transformation 



-1) 



k.Q: 



A2 



where 



?k 



,^.)V4+(A/3)2 



and 7k,p„ and r/k,p„ 
equation for A is 

1 



are quasi-particle operator*!^. The 



E 



tanh(e£°/3/2) 

Sk 



(22) 



the solution of which yields non-zero value of A and the 
critical temperature for having non-zero A. 

Next we shall examine whether there will be symmetry 
breaking of S also, for which we treat H' as perturbation 
on top of Hq, and the energy can be expanded in powers 
of S as Emf = Eo{A) + x(A)I]2 + 0(I]4). So whether 
there is an instability for developing finite S depends on 
the sign of x- Straightforward calculation shows that 



X(A) - - 



k,Q#/3 



2(e 



E 



(23) 

21^ 




rpp 



rpP 



■ N-M 



n 



N-H 



5 10 15 20 25 l^/En 

FIG. 2: (color on-line) Neel temperature Tjv as a function 
of optical lattice depth Vb. Four curves are p-band Tjv from 
Heisenberg model T^_jj, from mean-field theory T^_j^,j; s- 
band Tat from Heisenberg model T^_[j, from mean-field the- 
ory Tat-jv/, respectively. For parameters, we use '"'/S' in opti- 
cal lattices as an example^^. The lattice spacing A = 1064nm, 
K = n/X and Os — 15nm. The energy units is taken as Eb., 
which corresponds to 0.21fiK. 



where r,' 



a/3 



Pfi 



'^k.Pa '^k+TT.J 
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We numeri- 



cally compute Eq. (|23p and find that it is always positive, 
which indicates that there will be no symmetry breaking 
of E. Mathematically, it is because ■^k.pa ■Ck-i-Tf.p^ ex- 
cept for a few specific points in momentum space. In the 
limit of A ^ 0, the first term in the r.h.s. of Eq. 
diverges only around those a few points, while the sec- 
ond term diverges around the whole Fermi surface, hence 
the second term will always be dominative and the whole 
summation is positive. Physically, it is because the AF 
instability crucially relies on the nesting of Fermi surface, 
and the fact that there is no nesting between two differ- 
ent p-bands prevents inter-orbital spin order. Hence we 
have reached the second main result that 

In the weak interaction regime, the only spin order at low- 



(21) temperature is A = (Ei 



^ and this 



order coincides with that in strong interaction regime. 



V. NEEL TEMPERATURE 

Here we first calculate the wannicr functions in an op- 
tical lattice potential and use them to deduce the param- 
eters t, t' , U and W in our model. Then, we calculate the 
Neel temperature given by Heisenberg model of Eq. ([72]) 



Ttv-h, which is 1.276 Jox<5'(S' -I- 1)^ and also calculate 
Tn-m given by mean-field theory by solving equation 
((22)1 for the onset of non-zero A. We display our cal- 
culation of these two T/v for both half-filled p-band and 
s-band in Fig. O Tn-h decreases as the increase of the 
lattice height while T^-m increases as the increase of the 
lattice height. 

For small Vb, Tn-m 7n-h, the AF order is caused by 
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mean- field effect, while for large Vq where Tn-h *C 7n-m, 
the AF order is driven by super-exchange processes and 
the system is in strongly correlated regime. The regime 
where Tn-m and Tn-h cross is the crossover regime from 
weak to strong interactions, where both two methods 
will fail, and more advanced numerical technique like 
DMFT and QMC is required in order to give accurate 
Neel temperature^. We will leave this for further inves- 
tigation. However, the key point is that the s-band Tn-h 
is already of the order of O.lnK when the lattice height 
enters the strongly correlated regime of great experimen- 
tal interests, and it decreases to below 10~^nK rapidly as 
Vo increases, while the p-band Tn-h remains the order of 
InK even when Vq increases to ^ 20£'r. It is simply be- 
cause the wannier wave function in p-band is much more 
extended, its hopping amplitude t is one to two orders 
larger than that of s-band, which consequently leads to 
much larger Jox. 



In additional to higher Neel temperature, the spin-3/2 
of p-band AF state could have a better signal-to-noise 
radio than s-band spin-1/2 AF. All these provide great 
advantage for experimental study of quantum AF and 
its related phenomena due to super-exchange. Besides, 
we leave the effect of trapping potential to future inves- 
tigation which may introduce shell structure in density 
profile and the finite size effect into the problem. 

We hope that this work will open up a new route for 
studying strongly correlated Fermi gases in optical lat- 
tices, and be bases for future efforts along this direction. 
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